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ABSTRACT 

We present a study of the Milky Way disk and halo magnetic field, determined from observations of Faraday 
rotation measure (RM) towards 641 polarized extragalactic radio sources in the Galactic longitude range 100°- 
117°, within 30° of the Galactic plane. For \b\ < 15°, we observe a symmetric RM distribution about the 
Galactic plane. This is consistent with a disk field in the Perseus arm of even parity across the Galactic mid- 
plane. In the range 15°<|fc| < 30°, we find median rotation measures of — 15±4 rad m" 2 and -62±5 rad m" 2 
in the northern and southern Galactic hemispheres, respectively. If the RM distribution is a signature of the 
large-scale field parallel to the Galactic plane, this suggests that the halo magnetic field toward the outer Galaxy 
does not reverse direction across the mid-plane. The variation of RM as a function of Galactic latitude in this 
longitude range is such that RMs become more negative at larger \b\. This is consistent with an azimuthal 
magnetic field of strength 2 fj,G (7/^G) at a height 0.8-2 kpc above (below) the Galactic plane between the 
local and the Perseus spiral arm. We propose that the Milky Way could possess spiral-like halo magnetic fields 
similar to those observed in M5 1 . 

Subject headings: magnetic fields — Faraday rotation — polarization — Galaxy: halo 



1. INTRODUCTION 

Magnetism is an important component of the interstellar 
medium (ISM): it contributes to gas pressure to balance it 
against gravity, achieving hydrostatic equilibrium; it helps 
to remove angular momentum from a collapsing molecular 
cloud, allowing star formation to proceed; and it also confines 
and deflects cosmic rays. Yet, we still do not understand the 
details of processes that generate large-scale coherent galactic 
magnetic fields, especially those in regions outside the galac- 
tic disk. 

The existence of large-scale magnetic fields in galactic 
disks can be explained by the classical mean-field dynamo 
- the amplification of magnetic fields by small-scale turbu- 
lent motion (the a-effect) and di fferential rotatio n (the ui- 
effect) on a time scale of 10 9 years ( |Shukurov|2007| . Accord- 
ing to dynamo theory, an axisymmetnc magnetic field with 
quadrupole symmetry (even parity) with respect to the galac- 
tic mid-plane is most readily excited in the disk. This is found 
to be the case for the Milky Wa y - its disk field has even parity 
across the Galactic plane (e.g., Frick et al.|2001| Taylor et al. 
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2009; [Kronberg & Newton-McGee|201 T) . 

Dynamo modes excited in the thick disks or halos of galax- 
ies can have dipolar symmetry (odd parity) with respect to the 
galactic disk. These modes could also be oscillato ry, resulting 
in large-scale transient magnetic field reversals dSokoloff &| 



Shukurov 1990; Brandenburg et al.lr992[[Moss _et al .pHTUf . 
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Unfortunately, it is difficult to establish the halo field sym- 
metry with respect to the galactic plan e even in the most 
well-studied edge-on galaxy, NGC 253. |Braun et al.] pTTO) 
found that a large-scale magnetic field model of an axisym- 
metric spiral form plus an out-of -plane quadrupole extension 
that could reproduce polarization patterns observed in a num- 
ber of face-on galaxies. Unfortunately, the backside halos in 
their target galaxies were severely depolarized and the field 
symmetry with respect to the mid-plane could not be deter- 
mined. Disk and halo fields could be of different axial ge- 
ometries: for example, M 51 (inclination of -20°) was found 
to ha ve a bisymmetric ha lo field and an axisymmetric disk 
field dFletcher et al.|2011| >. 

Characterizing the halo magnetic field symmetry in our own 
Milky Way is even more challenging because of our location 
within it. Early works based on a sparsely sampled all-sky 
Faraday rotation measure (RM) distribution of pulsars and 
extragalactic sources (EGSs) have suggested a dipolar-type 
Galactic halo ma gnetic field d Andreas sian & Makarov 1988 



|Han et al.||1997| l. However, using the RM of ~ 1,000 ex 
tragalactic sources, |Mao et aT7| d2010[ l demonstrated that the 
Milky Way's vertical magnetic field at the location of the 
Sun is of neit her quadrupolar n or dipolar symmetry. More 
recently, Wolleben et al. ( 2010| l used RMs of diffuse polar- 
ized Galactic synchrotron emission to show that the EGS RM 
pattern thought to be the signature of a dipolar halo field is 
at least partly due to a local magnetized shell. Clearly, fur- 
ther effort is needed to properly describe the geometry of the 
Galactic halo magnetic field. 

In this paper, we characterize the symmetry of the Galac- 
tic halo magnetic field by studying RMs of EGSs in the 
region 100°<Z<117°, |b|<30° (towards the Perseus Arm). 
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This study is based on existing EGS RMs in the Cana- 
dian Galactic Plane Survey (CGPS) latitude extension region 
(100°</<117°, -3°<Z?<+18°, Brown et al. in prep) and a 
new Very Large Array (VLA) RM survey that extends the lat- 
itude coverage of the RM data to \b\ <30°. The latitude cover- 
age of this new VLA data set enables one to study the parity of 
the halo field across the Galactic plane. The longitude range is 
ideal because it samples the magneto-ionic medium half way 
between the two strips of the Galactic plane that will be cov- 
ered by the GALFA C ontinuum Transit Survey (GALFACTS) 
flaylor & Salter|2010| i. 

In Section |2J we describe the VLA observations and data 
reduction. We then compute EGS RMs and present their spa- 
tial distributions in Section [3] In Section HI we infer proper- 
ties of the symmetric disk field towards thePerseus spiral arm 
using the observed R M-v s-fc trend wi thin ~T5° of the Galac- 
tic plane. In Section |5TTj and Section [572] we show that RMs 
in the surveyed region are likely produced by a large-scale 
ma gneti c field in the warm ionized medium (WIM). In Sec- 
tion 5.3 we illustrate that the existing halo field models can- 



not reproduce satisfactorily the RM-vs-Z? trend in our survey 
area . We construct a simple Galactic halo field model in Sec- 
tion [5j4] that can successfully reproduce the general observed 
RM trend in the surveyed area. 

2. OBSERVATIONS AND DATA REDUCTION 

To directly measure the parity of the Milky Way's disk and 
halo magnetic field in a latitude strip, we extend the latitude 
coverage of the CGPS latitude extension region by conduct- 
ing an EGS RM survey with the VLA. We construct a mag- 
netic field model for the Milky Way to estimate the number 
of sources needed for such a detection. We assume the Milky 
Way has a symmetric 1 /iG disk field with a scale height of 1 
kpc and an antisymmetric 1 /iG halo field with a scale height 
of 4 kpc. In additio n, we use the therma l electron thick disk 
model derived by Corde s~& Lazio| ( |2002| l as the Faraday ro- 
tating layer, for the purpose of observation planning. We find 
at latitudes higher than 30°, the difference in RMs above and 
below the plane due to the halo field might be too small to de- 
termine its parity, hence we use |Z?|=30° as the latitude bound- 
ary of our observations. This halo field toy model predicts 
an RM difference of ~ 40 rad m" 2 for > 10°. Since the 
standard deviation of EGS RMs at |fc|~30° is ~ 18 rad m" 2 



(e.g.,|Schnitzeler 2010), one needs approximately 15 sources 
in each 1° latitude bin to distinguish the RM trend for b>0° 
from b<0°, resulting in a source requirement of ~ 900 EGSs 
within 30° of the Galactic plane. The CGPS latitude exten- 
sion survey provides ~ 300 extragalactic sources in the lati- 
tude range -3.5°<£><+18° (Brown et al. in prep). Therefore, 
we have selected ~ 600 compact extragalactic sources from 
the N RAO VLA Sky Survey (NVSS) catalogue ( jCondon et al. 
1998) with polarized intensity brighter than 4 mJy in the re- 
gion 100°<Z<117°, -30°<£<-3° and +18°<Z?<+30° as tar- 
gets. 

Observations of the selected polarized EGSs were con- 
ducted at the VLA on 2008 June 8th and 23rd, July 7th and 
8th, August 26th and 31st, and September 1st and 2nd in D 
and DnC configurations. Data were taken in two 25 MHz 
bands in the spectral line mode. Each band consists of seven 
3.125 MHz-wide channels (with an edge channel discarded), 
centered at 1365 MHz and 1486 MHz, respectively. Data 
taken on June 8th had the higher frequency band centered at 
1512 MHz. Unfortunately, these latter data were rendered un- 
usable due to severe radio frequency interference (RFI). The 



two frequency bands are separated by 100 MHz to ensure 
precise RM determination, as the uncertainty in RM is in- 
versely proportional to the span in the coverage in wavelength 
squared. Each source was observed for ^1.6 minutes at each 
frequency band. Including slew times and overheads, the to- 
tal observing time was ~ 40 hours. Standard VLA primary 
calibrators 3C147 and 3C286 were observed at the beginning 
and the end of each observing run for absolute flux and po- 
larization leakage calibration. Polarized sources 3C286 and 
3C48 were observed for the purpose of absolute polarization 
angle calibration. In addition, a secondary calibrator (one of 
the following: 1800+784, 2005+778, 2133+826, 2355+498, 
0029+349) located closest to the target EGSs was observed 
every hour to obtain time-dependent antenna gains. 

We converted the raw data into UVFITS format using the 
task FITTP in AIPS in order to carry out data calibration, 
and then m ake images using the MIRIAD package ( [Sault l 

id 



et al. 



1995 Flux densities of the sources were calibrated 



to 3C147, whose Stokes I value at 1365 MHz and 1486 MHz 
is 22.5 Jy and 21.2 Jy, respectively. The absolute polarization 
angles of the sources were calibrated to either 3C286, whose 
angle at 1365 MHz and 1486 MHz is +33°, measured from 
North through East, with a rotation measur^Jof rad m" 2 or 
3C48, whose angle is tied to that of 3C286: -55° at 1365 
MHz and -30° at 1486 MHz. Time-dependent antenna gains, 
polarization leakages and absolute phases of the 14 spectral 
channels were calibrated individually. For each EGS, images 
of Stokes Q and U were made in each of the 14 3.125-MHz 
wide frequency channel using natural weighting to maximize 
sensitivity. This results in channels maps each with a sen- 
sitivity of roughly 0.4 mJy/beam. We excluded the shortest 
baselines (uv distance < 0.5 kA) to avoid imaging extended 
Galactic synchrotron emission. The Stokes Q and U channel 
maps were then deconvolved and restored to the resolution at 
1356 MHz (~1'). A linearly polarized intensity (PI) map cor- 
rected for positive bias was made for each source. The bright- 
est polarized pixel of each EGS was identified and its Stokes 
Q and U values across the frequency band were extracted for 
RM determination. 

This work also makes use of 339 EGS RMs in the CGPS 
latitude extension region (100°<Z<1 17°, -3.5°<£<+18°, 
Brown et al. in prep). These EGSs have a mean total in- 
tensity of 81 mJy/beam and a mean polarized fraction of 5%. 
RMs are determined using Stokes Q and U measured in four 
7.5-MHz channels centered at 1420 MHz. The typical RM 
error is 15 rad m" 2 . 

3. RESULTS 

3.1. RM Computation 

Faraday rotation is a birefringence effect (see e.g. Section 
4.3 of |Mao et al.||2010| for more details). We have com- 
puted rotation measures using RM Synthesis and R MCLEAN 

~ 2009] l, following 
). The RM uncer- 



( |Brentjens & de Br uyn 200 5] |Heald et al 
the algorithm presented by |Mao et al. | (j20 1 



tainty is computed by directly measuring the noise in the real 
and imaginary parts of the RM spectrum. Our VLA observing 
frequency setup results in RM spectra with FWHM ~ 400 rad 
2 and a maximum detectable RM of 10 4 rad m~ 2 . 



m 



Recently, |Farnsworth et al. ( 201 l| l demonstrated that using 
RM synthesis alone might not be sufficient to determine the 

1 See also http://www.ph.unimelb.edu.au/~rsault/software/miriad/vla- 
polarimetry.html tor specifics of reducing VLA data in MIRIAD. 
z http://www.vla.nrao.edu/astro/calib/manual/ 
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underlying Faraday structure, even in the simple case of two 
RM components with different intrinsic polarization angles. 
This is because RM Synthesis does not have an equivalent 
of a reduced x 2 to measure the goodness of fit as in the least 
square fit of polarization position angle as a function of A 2 . As 
a result, the solution can converge to an incorrect RM value. 
To ensure the reliability of our RMs, we have computed the 
reduced \ 2 of the angle against A 2 relation using the RM value 
obtained from RM synthesis, and we only accept RMs with \ 2 
< 2. In addition, we have inspected the behavior of Stokes Q, 
U and the polarized intensity as a function of A 2 for each EGS 
to ensure that the result from RM synthesis and least-square 
fit are in reasonable agreement. We have derived reliable RMs 
for 302 EGSs in the surveyed region, as listed in Table [T] 

3.2. Comparison with NVSS RMs 

[Taylor et aL] ( [2009] ) have computed RMs of EGSs in the 
NVSS catalogue using the data from the original surve y. We 
find matches fo r 244 sources of our VLA RMs in the |Tay-| 
|lor et al.| ( |2009[ l catalogue. In Figu re [T] we have plotted the 
RMs derived by |Taylor et al.| ( |2009| l against the RMs that we 
have derived for the same sources. The solid line of slope 1 
indicates where sources should lie if the NVSS RMs and our 
RMs are equal. Approximately 57% of the RMs from the two 
samples agree with each other within their measurement er- 
rors. The linear correlation coefficient between the two RM 
data sets is 0.47. However, if we compute the coefficient after 
discarding 3 EGS with the most extreme |RM| in the matched 
data set, the linear correlation coefficient increases to 0.89, 
suggesting good RM agreements between the two data sets. 

3.3. RM Distribution in the Surveyed Region 

The distribution of EGS RMs in our VLA survey and in the 
CGPS RM survey (Brown et al. 2003, Brown et al. in prep) in 
the Galactic longitude r ange 100 o -117° an d latitude |b|<30° 
are over-plotted on the [Pinkbeiner ( 2003 i all sky Ha com- 
posite map in Figure [2] RMs are large and negative for sight 
lines within 10° of the Galactic plane, then |RM| slowly de- 
creases with Galactic latitude to |Z?|^10°. At higher latitudes 
(|£>|>15°), the magnitude of RM toward the northern Galactic 
hemisphere remains small but negative, while RMs toward the 
southern Galactic hemisphere slowly become more negative. 

This RM against b trend is illustrated in Figure [3] In the 
top panel, median RMs are plotted as a function of b us- 
ing latitude bins containing equal number of EGS each, with 
error bars denoting standard error of the mean within each 
bin. To highlight the difference between the behavior of RM 
above and below the plane at \b\ < 15°, we have folded the top 
panel about b=Q°, producing the figure in the bottom panel. 
This plot demonstrates that the RMs are symmetrically dis- 
tributed close to the Galactic plane. The difference in RM for 
sight lines above and below the plane becomes pronounced at 
|b|>15°. 

4. THE SYMMETRIC LARGE-SCALE MAGNETIC FIELD IN THE 
DISK TOWARDS THE PERSEUS ARM 

Magneto-ionic structures along the sight lines toward the 
outer Galaxy are less complex than those towards the inner 
Galaxy as the former only pass through one major spiral arm: 
the Perseus arm. For this reason and because of the high 
angular density of RM measurements in our new sample at 
100°<Z<117°, the symmetry properties of the disk field can 
be revealed by plotting RM against b (Figure|3| without need- 
ing to apply wavelet analysis (Frick et al. 20u~l ) or averaging 



algorithms (Kronberg & Newton-McGee 2011 ). The signs of 
RM above and below the disk within ~ 15° of the Galactic 
plane are the same, implying that towards the outer Galaxy, 
the plane-parallel disk magnetic field preserves its direction 
across the Galactic mid-plane (even parity), in agreement with 
previous works. This picture fits in well with our current un- 
derstanding of the disk field origin: the dynamo mechanism 
predicts a symmetri c disk field as it is the ea siest mode to be 
e xcited (see e.g.,|Wid row 2002; Shukurov 2007). 

|Rae & Brown|( |2010[ ) investigated the symmetry plane of 
the disk field by fitting a Gaussian to the RM versus b be- 
havior using EGS RMs in the CGPS latitude extension region 
which mostly probe sight lines above the plane (latitude cov- 
erage -3.5°<b<+18°). They found that the RM distribution 
is not symmetric about b=0° , but is rather shifted to the north 
by 1°. Rae & Brown (2010) attributed this offset to a warped 
Galactic magneto-ionic disk. A warp in the Galactic disk to- 
wards the outer Galaxy has b een detecte d in the stellar com- 
ponent (for a su mmary, see |Vallee|20 11) as well as in neutral 
hydrogen ( Kalberla & Kerp 2009 ). The distance to the warp is 
estimated to be at least 5 kpc from the Sun towards the outer 
Galaxy (e.g., |Levine et al.|2Q08) >. Our new VLA observations 
provide EGS RMs towards negative latitudes in the same lon- 
gitude range as the CGP S data. This pu ts us at a position to 
re-examine results of|Rae & Brown| ( |2T)T0| ). 

Similar to Rae & Brown| ( |2010| l, we center the RM distribu- 
tion about the Galactic disk by fitting a Gaussian^] The Gaus- 
sian centroid fit is sensitive to the latitude range within which 
RMs are fitted. We use |Z?|=10° as the initial boundary of the 
fit since it is where the RM-vs-Z? trend changes abruptly. This 
is likely because RMs of EGSs at higher Galactic latitudes are 
dominated by the Galactic halo magnetic field which has a dif- 
ferent vertical symmetry from the disk field. Including RMs 
at higher b when fitting will then likely degrade the quality 
of the centroid fit. We vary the latitude boundary of the fit 
until the reduced \ 2 of the Gaussian fit is minimized. We 
found that the distribution of rotation measure is symmetric 
about /?=-0.13°±0.48°, with corresponding fit boundaries at 
£=±22.15° . This is different from the offset of the RM sym- 
metry axis to positive b reported in Rae & Brown ( 2010). This 
is likely due to the limited Galactic latitude coverage (mostly 
above the Galactic plane) of the CGPS data, and the fact that 
we have used a rece ntly revised CGPS R M catalogue, differ- 
ent from that used in Rae & Brown (2010). Our findings of the 
RM symmetry with respect to the Galactic mid-plane at low 
G alactic latitude is in agre ement with the conclusion reached 
by Lazio & Cordes (1998), who showed, using the scattering 
of EGSs, that the ionized Galactic disk is not warped towards 
the outer Galaxy. 

The fitted RM distribution has a width of 5. 49° ±0.78°: the 
|RM| produced by the disk field is less than 5 rad m" 2 for 
\b\ >15°. The symmetric RM behavior about the Galactic 
plane within ~ 15° is con sistent with the result of |Kronberg| 
|& Newton-McGee] ( |201 l| l. These authors used a smoothed 
all sky RM data set to show that the Milky Way's disk mag- 
netic field stays symmetric about the mid-plane for |£>|<15°. 
We note that translating this latitude boundary into a physi- 
cal height below which the disk field dominates requires the 
knowledge of where along the line of sight most of the Fara- 
day rotation occurs. If we assume that most Faraday rotation 

3 A parabolic fit to the RM distribution provides a centroid location con- 
sistent with that from a gaussian fit. However, we chose to use the Gaussian 
fit since it gives a lower reduced \ 2 . 
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is mainly produced by magnetic fields and electrons in th e 
Perseus spiral arm at a distance of 2 kpc ( Xu et aLp 006b), 
then the disk magnetic field dominates within ~ 540 pc of the 
Galactic plane. 

5. THE ROTATION MEASURE DISTRIBUTION TOWARDS 

100 o <L<117°, 15°<|B|<30° 

In this section, we investigate the RM behavior in our sur- 
veyed longitude range at 15°<|b|<30°. We first establish the 
location of the Faraday rotating layer and then determine if 
structures in thermal electrons or in magnetic fields are re- 
sponsible for producing the observed RMs. Finally, we con- 
struct a simple model of the Galactic halo magnetic field to 
explain the observed RM trend as a function of Galactic lati- 
tude. 

5.1. The location of the Faraday rotating layer at |Z?|>75° 

Since Faraday rotation is an integral from the EGS to the 
observer, the observed RM could originate anywhere along 
the line of sight. Therefore, it is important to determine where 
most of the Faraday rotation takes place. The observed RM 
could be of extragalactic or Galactic origin. Galaxy clusters 
and superclusters are magnetized and can produce enha nced 
RMs in directions away from the Galactic plane (s ee e.g., |Kim| 
|eTaLl[T990l fClarke et al.|[200T] |Xu et aLl|2006a| i. SincelEe 
mid-Galactic latitude region in our VLA survey does not co- 
incide with extended galaxy clusters or superclusters, we rule 
out the possibility that the RMs we measure toward the survey 
area are produced in the intracluster medium. As h igh veloc- 
ity clouds (HVCs) c an produce observable RMs (McClure- 
|Griffiths et al.|20i0) l, we verified, by co mparing maps of the 



high velocity neutral hydrogen sky (e.g., Wakker & van Wo- 
erden|1997] ), that no HVC lies within our surveyed region. 



Flie negative RMs seen towards b^-2Q° in the longitude 
range 100°- 117° in Figure [2] belong to Region A, an anoma- 
lous RM region (60 o <Z<140°, -40 o <fr<+10°) ([Michel & | 
Yahil|19731 Simard-Normandi n & Kron berg 198QJ>. JSimard^l 
Normandin & Kronberg| (1980) suggested that the Sun does 
not reside in the region producing the negative RMs due to 
the lack of a positive RM region in the opposite direction. 
Since Region A appears to be located wi thin the boundary of 
Galactic radio Loop II (the Cetus Arc), Berkhu ijien] ( ] 1 97 i) 
suggested that the unusually large and negative RM in Region 
A is produced by Loop II. Loop II is thought to be the shell 
of a nearby (~ 100 pc) supernova remnan t expanding into 
the ma gnetized interstellar medium, as by Spoelstra ( |1972| l 
and by |Vallee| ( [T982l [19931 . m their models, a shell of thlck- 
ness 10 pc, with an electron density of 1.2 cm" 3 and a mag- 
netic field strength of ~ 10 could produce RMs of simi- 
lar magnitudes to those observed within Region A. However, 
one expects the thermal electrons in this shell to produce an 
emission measure of ~ 14 pc cm" 6 . This is inconsistent with 
WHAM observations: the mean emission measure in Region 
A is roughly 2 pc cm" 6 , far smaller than the predicted value 
(since region A is far outside the Galactic plane, we expect 
little extinction). Furthermore, there is no enhancement in 
Ha along the boundary of Region A, that might correspond 
to limb-brightening. One alternative is that a magneto-ionic 
medium which is non-Ha emitting i s responsible for pro- 
ducing the observed Faraday rotation. Simard-Normandin & 



cosmic rays was too low to explain the observ ed RM mag- 
nitude. Simard-Normandin & Kronberg (1980) further sug- 
gested that given the dimension of Loop II, the total energy 
injection requires multiple supernova explosions. It is unclear 
how magnetic fields can stay coherent on such a large area in 
the sky after a few of these violent events unless they produce 
a single coherent superbubble. We also note that the boundary 
of Loop II does not lie within our survey area. 

More re cently, |Stil et al.| (|2011[ ) used RMs derived from 
the NVSS (Tay lor et al.|2009| l to revisit the possible connec- 
tion between Loop II and Region A. They argued that the 
RM sign-change from negative within to positive outside Re- 
gion A appears to coincide with the bound ary of Loop II at 
/^160°. This was taken as evidence by |Stil et al. (2011 1 of 
a magnetic field reversal due to Loop II. Upon close inspec- 
tion, however, we find that positive RMs outside Region A 
near Z^160° is likely produced by ionized gas and magnetic 
fields in the HII region Sh2-220 (The California Nebula) and a 
more extended Ha e mission region around it (Reynold s|1988[ 
Haffn eret al.|| 1999 1. The coherent magnetic field prop erties 
of Sh2-220 have been investigated in detail by Harvey-Smith| 
et al.| ( |20TT] >. We suggest that the RM sign change is likely at 
least partly due to the presence of the HII region, rather than 
due to the absence of Loop II. 



Finally, using RMs towards pulsars in Region A, Newton 



Kronberg (1980) considered the possibility that low energy 



McGee (2009i reached the conclusion that Faraday rotation 
along sight lines within Region A must be produced at least 
0.95 kpc away from the Sun. This once again weakens the 
possible connection between Region A and Loop II, the latter 
of which is thought to be located within ^100 pc. There- 
fore, we suggest that the negative RMs that we see at mid- 
latitude below the Galactic plane are likely not produced by 
local structures in the ISM, but rather by larger scale struc- 
tures located at least 1 kpc from the Sun. 

Our surveyed region above the Galactic plane lies within 
the boundary of Radio Loop III (Berkhuijsen 1971 ). Although 
RMs at the low Galactic longitude edge of Loop III appear to 
be enhanced (see Figure [4), the RM distribution within Loop 
III itself does not appear to be anomalou^] 

5.2. Enhancement in Electron Density or Magnetic Field? 

Since RM is sensitive to magnetic field weighted by ther- 
mal electron density, an increase in EGS RM could reflect 
either an increase in the average line-of-sight magnetic field 
strength, in the thermal electron density or in both quanti- 
ties. If thermal electrons are responsible for pro ducing the 
observed RM as Simard-Normandin & Kronberg] ( |1980| > ar- 
gued, and the same thermal electrons also emit in Ha, then 
one can distinguish between RM due to an increase in line- 
of-sight magnetic field strength from that due to an increase 
in electron density by inspecting the distribution of Ha emis- 
sion along the sight lines of interest. 

We found that for both local (-25 km s" 1 <Vlss<+25 
km s" 1 ) and Perseus arm (-60 km s~ 1 <Vlsr < ~40 km 
s" 1 ) ionized gas, the emission measure is distributed sym- 
metrically about the Galactic plane - there is comparable 
amount of thermal electrons at +20°<b<+30° as that to- 
wards -30°<£><-20°. Therefore, the structure in RM is un- 
likely to be due to enhancements in thermal electrons. In 
stead, it is likely that the observed RM in the surve yed re 
gion reflects changes in magnetic field structures. Simard 



relativistic electrons, which are not visible in Ha, produce the 
observed RM, but the estimated density of these low-energy 



4 We note that the variance_^fRM_^Wun_Ij00^riI^as^e£orte^^o^e 
enhanced by Stil et al.MZUlll. 
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Normandin & Kronberg ( 1980[ > also pointed out that if the 
increase in RM was primarily due to an increased electron 
density, then the electron density would be a factor of 50 
larger towards the southern Galactic hemisphere than towards 
the northern Galactic hemisphere (this is roughly the ratio of 
|RM| at b=-20° to that at b=+20°). If the pulsars were be- 
yond the Faraday rotating region, their DMs would reflect this 
increase in electron densities. If that is the case, then pul- 
sars at b>+20° should on a verage have much smaller DMs 
than those at b<-20°. Both Simard-Normandin & Kronberg 
(1980) and |Newton-McGee| ( f2009| ) did not find evidence for 
such a pulsar DM distribution, supporting our claim that RMs 
in our VLA survey likely reflect structures in magnetic fields 
rather than electron densities. 

5.3. Prediction of RM from Existing Halo Magnetic Field 
Models 



In Sections |5.1| and |5.2| we have established that RMs to- 
wards mid-latitudes in the longitude range 100°- 117° are of 
Galactic origin and are likely due to large-scale structures. 
We also reach the conclusion that RMs at mid-Galactic lat- 
itudes likely reflect structures in magnetic fields rather than 
electron densities. In this section, using expressions for the 
halo magnetic field in the literature, we show that no existing 
models are successful in reproducing the observed EGS RMs 
for sight lines |Z?|>15° in the longitude range 100°<Z<117°. 
It is evident that additional modeling is required to properly 
characterize the Galactic halo magnetic field. 

The median RM for b>+\5° is -15±4 rad m" 2 while the 
median RM for b<— 15° is -62±5 rad m" 2 . Since the me- 
dian RM is negative both above and below the plane at mid- 
Galactic latitudes (|Z?|>15°), we suggest that the large-scale 
halo magnetic field parallel to the plane toward the outer 
Galaxy does not reverse direction across the mid-plane. This 
is consistent with recent findings of Pav el et aLl p012 ): these 
authors used Near IR polarimetry of starlight towards / = 150° 
between -75° <b< +10° and simulations to rule out mag- 
netic fields of odd parity towards the outer Galaxy. Since we 
did not find an asymmetric electron density di strib ution about 
the Galactic mid-plane at |b|>20° in Section |5.2| the differ- 



ent RM magnitudes or orientations above and below the plane 
are likely due to different halo magnetic field strengths in the 
northern/southern hemispheres. 

Expressions for the halo magnetic field in the literature of- 
ten contain a magnetic field reversal across the Galactic plane. 
|Han et al.| (1997} and |Han| P002| ) claim, based on the anti- 
symmetric sign of EGS RMs at high Galactic latitudes to- 
wards the inner Galaxy, that the Milky Way's halo magnetic 
field is dipole-like (for which the horizontal component re- 
verses direction across mid-plane, but t he vertical component 
does not). This claim is challenged by |Mao et al. ( 2010j ), who 
demonstrated that the magnetic field geometry at the location 
of the Sun canno t be a pure dipole or a quadrupole. Moreover, 
|Wolleben et al.| ( |2010| ) used RMs determined from diffuse 
synchrotron polarized emission to show that a nearby mag- 
netized bubble could produce p art of the anti-symmetric EGS 
RM pattern that Han et al.] ( |1997| ) interpreted as a signature of 
a global dipole-like field. Although models in the literature 
that predict a reversal in magnetic field direction across the 
plane towards Galactic quadrant 2 are unlikely to be physi- 
cal in light these more recent works, we nonetheless plot the 
predictions from these models in this section. 

Existing halo magnetic field m odels consist of two gen- 
eral types: toroidal fields (e.g., |Sun et al.|200"8] ) or exponen- 



tial plane-parallel halo fields that fall off with Galacto-centric 
radius and distance from the mid-plane. Parameters of the 
halo field models can be divided into two different groups as 
well: values that can qualitatively reproduce the halo field ob- 
servables (RMs of EGSs, diffuse synch rotron emission and 
starlight polarization)(e.g., |Sun et al.|2008[ ), or values that are 
obtained by finding a best fit to the halo fi eld observables, 
minimizing the \ 2 ( e -E-- 



2009). 



, Jansson et al. 

In the following, we list all the halo field expressions in the 
literature. We define a parameter S that captures the dipolar 
nature of the halo field: its value is -1 for b>0° and +1 for 
b<0°. The sign is chosen such that the halo field points in the 
same direction (clockwise) as the disk field below the plane 
and reverses such that it is (counterclockwise) above the plane 
in Galactic quadrant 2 when viewed from the north Galac- 
tic pole. The halo field expressions are given in a cylindrical 
Galacto-centric coordinate system, for which components in 
the radial (r), azimuth ((j>) and vertical (z) directions are listed 
separately. We also provide the corresponding model parame- 
ters and note whether these parameters are found qualitatively 
by ey e or by a qu antitative fit to halo field observables. 
The Jansson et al. (2009 ) halo field (£ , 2mi,2009) is given by: 



For r < r c 

B r = SB"sin(p inner ), 
B 4> = -SB"cos(p inner ), 

By = 0. 

For r > r c 

B r = B"sin(p ouler ), 
B c f > = -B"cos(p outer ), 
B z = 0. 



(1) 



The best fit parameters are: r e =8.72 kpc, the radius beyond 
which the horizontal field reverses across the mid-plane; Bq 
= 2.3 fiG, the strength of the horizontal field; and p mner =-2° 
(Pouter = -30°), the pi tch angle for r < r^ (r > r c ). 
A variation of the Prouza & Smfda (2003) double torus 



magnetic field of the following form has been heavily used 
in the literature: 



B,- = 0, 

B<t> = SBq - 

B z = 0. 



1 + 



1 



2 r H 



expi- 



re 



(2) 



The values of B%, 



4>Z 



and zf used by 



Sun et al. 



(2008) 



Jansson et al. (2009), Sun & Reich (2010) and Pshirkov et al. 
( 20 1 1 [ ) are summarized in Table |2| The nature of these pa- 
rameters (qualitative or fit) are listed in the last column of 
Table [2] To illustrate the magnetic field strength of the dou- 
ble torus model, we have plotted the azi muthal magnet i c field 
strength using the best fit parameters of Jansson et al. ( 2009 ) 
in a vertical cross-section of the Milky Way in the top panel 

of Figure [5] 

Finally, "R uiz-Granados et aL] ( |2010| ) used a bi-toroidal halo 
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field: 



B, = 0; 

Ba, = — l - arctan(—)exp(- 

rl+r (X\ 

B- = 0.2. 



2a\' 



(3) 

The corresponding azimuthal magnetic field strength in a ver- 
tical cross-section of the Milky Way is shown in the bottom 
panel of Figure [5] The best fit values are r\ > 33.8 kpc, <j\ = 
2.9 kpc and > 4.7 kpc. 

In computing expected RMs from various halo field mod- 
els, we assume t hat the Sun is l ocated R slm = 8.4 kpc from the 
Galactic center (Rei d et al.|2009| l. We further assume a disk- 
field-to-halo-field transition heighj^Jof ~ 540pc (as inferred 
from the RM-vs-Z? trend presented in Section|4]i. For the Fara- 
day rotating electrons, we use the thermal electron thick disk 
model of |Gaensler et al.| ( [2008] l, in which the WIM is mod- 
eled as an exponential disk of scale height /zwim=1-83 kpc 
with a mid-plane density near the Sun of n e .o=0.014 cm" 3 . We 



adopt the radial dependence of the thick disk used by Cordes 
|& Lazio| ([2002 ) with a truncation at Galacto-centric radius 
Ai= 20 kpc. The exact electron density distribution used is 



n e = n e n 



exp(- 



(4) 



Cordes & Lazio (2002i adopt a vertical distribution in the 
form of sech z (z/hwiM) rather than the exponential in the above 
expression. We verify that the two z dependence produce RM 
against Galactic latitude trends that are very similar. Since 
we are focusing on fitting mid-latitude RMs, we do not seek 
to fit to EGS RMs at |fe|<15°. We numerically integrate 
the thermal electron density weighted magnetic fields toward 
Z=108.5°, the median longitude of our surveyed area. The pre- 
dicted RM-vs-Z? trends are shown in Figure [6] We note that a 
low disk-to-halo transition height is strongly favored simply 
because there are more thermal electrons closer to the mid- 
plane, which produce larger RMs given the same magnetic 
field distribution. If one adopts larger transition heights, the 
predicted |RM| will be smaller than that shown in Figure [6] 
due to the decrease in thermal electron density at larger dis- 
tances from the Galactic plane. 

It is clear from Figure[6]that existing Milky Way halo mag- 
netic field models do poorly in reproducing the observed RM- 
\s-b trend for both positive and negative latitudes in the sur- 
veyed region. None of the models can simultaneously re- 
produce the symmetry property of the observed RM and the 
magnitude of the RMs. The disagreement between observed 
RMs and model predictions i n this / and b r ange hav e been 
noted by sev eral authors. Both |Sun et al.| ( [2008] ) and Pshirkov 
et al. ( |201 \\ have attributed the discrepancy to radio contin 



uum loops and Region A, but we have shown in Section 5.1 
that these are unlikely to be the cause. 

5.4. A Model of the Milky Way Halo Field towards 
100°<l<117° 

We first attempt to vary parameters of the existing halo field 
models to match our observed RMs. Unfortunately, adjusting 

5 Below this height, the disk magnetic field dominates and the coherent 
halo field is assumed to be zero. Above this height, the halo magnetic field 
dominates and the coherent disk field is assumed to be zero. 



the parameters to match the observed RM magnitudes pro- 
duces excess RM towards the Galactic poles, for which the 
observed values are zero towards the north and +6 ra d m~ 2 to- 
wards the south flay lor et al.|2009||Mao et al.|201Q) . Instead 
of fitting the observed RM with a global exponential field or 
a double toroid above and below the plane, we choose to ap- 
proach the problem by constructing a very simple model of a 
part of the Milky Way's magnetic halo to reproduce the ob- 
served RM trend. 

Since we reside in the Milky Way, drawing conclusions on 
the Galactic halo magnetic field configuration based on EGS 
RMs is challenging without a bird's eye view. Observations 
of external galaxies ha ve revealed X-shap ed magnetic fields in 
edge-on galaxies (e.g., Heesen et al.|2009} and axisymmetric 
spiral with an out-of-plane quadrupole extension in face-on 
spiral galaxies (e.g., Br aun et al.p OlO). None of these ex- 
ternal galaxies possess magnetic field structure similar to the 
Milky Way halo field models proposed in the literature: a pure 
toroidal field or a simple exponential field. Therefore, current 
Galactic halo field models might be unrealistic. 

On the theoretical front, the popular double-torus model 
might be physically too simplistic as its pitch angle is zero 
everywhere in the halo. It is thought that a pure circular mag- 
netic field cannot be maintain ed by any veloc ity field against 
turbulent magnetic diffusion (Beck et al .|1996[ l. Furthermore, 
if the halo magnetic fields were of dynamo origin, then spiral 
magnetic field lines would form naturally, since the dynamo 
mechanism relies on differential rotation to transform radial 
fields into azimuthal fields. Motivated by both observation 
and theory, we seek a new halo magnetic field model to ex- 
plain the RMs in our observed region. 

We assume that there is a constant, coherent azimuthal 
magnetic field Bq parallel to the Galactic plane in Galacto- 
centric radius range r; nner < r < r outei kpc and at a height 
Ziower > \z\ > z U pper kpc from the Galactic mid-plane. We as- 
sume the halo magnetic field strength outside this specified 
region is zero: 

Foi t"[ nner <C T <C Pouter; 

Slower ^ |^| ^ Supper • 

B r = 0, 

Ba, = —Bk, 



B z = 0; 
otherwise : 

B r = B^=B z = 0. 



(5) 



We allow different values of Bq above and below the mid- 
plane and assu me that the Faraday rotating medium is the 
WIM thick disk ( |Gaensler et al.|2008) . 

We found by performing a search in parameter space that a 
clockwise toroidal halo field with By = 2 /iG above the plane 
and Bq = 7 /iG below the plane, rj nner = 8.8 kpc, r outer = 10.3 
kpc, Ziower=0.8 kpc and z upP er = 2 kpc can successfully repro- 
duce the RM-against-/? trend in our observed area for \b\ >15°. 
This is demonstrated in Figure U\ where the predicted RM 
variation as a function of latitude towards Z=108.5° (the av- 
erage Galactic longitude of the observed region) is plotted. 
The shape of the model-predicted RM closely matches that 
of the data. The 2D spatial distribution of the predicted RMs 
in the surveyed area is plotted in the left panel of Figure [8] 
where EGS RMs at |fe|>15° are overlaid as circles with the 
same color scale. The color of filled circles and that of the 
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background (prediction from our simple model) are in gen- 
eral agreement. The right panel of Figure j8] shows the distri- 
bution of the residual RMs after subtracting the model pre- 
diction from the observed RMs. The residual RMs are small 
and mostly consistent with zero. We have compared the re- 
duce x 2 of our simple model to all of the models presented 
in Section [53] at |i|>15°, in the Galactic longitude range 
100°</<1 17 u . Our simple model is superior to existing halo 
field models at larger than 99.9 % confidence level in the re- 
gion of interesfj 

The modeled coherent magnetic field exists at a distance 
between 0.8 and 2 kpc from the Galactic disk, which is well 
above the z height where the disk field dominates (~ 540 
pc, Section |4|. Hydrostatic equilibrium in the Galactic halo 
predicts a magnetic field strength of ~ 3-4 /xG at a dis- 
tance 1-2 kpc from the mid-plane fKalberla & Kerp]|1998[ 
IFletcher & Shukurov|2qoT] |Cox|2005) . The best-fit magnetic 
field strengths towards the northern and the southern Galactic 
hemispheres of 2 and 7 /iG are in rough agreement with hy- 
drostatic equilibrium. The need for a strong (7/iG) field could 
be a consequence of the simplicity of our model: we have 
set the coherent magnetic field strength to be zero external to 
Galactocentric radii in the range 8.8-10.3 kpc. Although our 
best-fit value for the magnetic field strength in the southern 
halo is comparable to that derived in the double-torus model 



used by Sun et al. ( 2008[ ), our model is unlikely to o ver-predict 
the cor responding synchrotron emission as did the |Sun et aL 
( |2008[ ) model. This is because our 7/xG halo ma gnetic field 
has a much lower volume filling factor than the |Sun et aL 
(2008) double-torus field model. In our model, the distance 
to the region where Faraday rotation occurs is approximately 
1 k pc, si milar to the limit implied from pulsar RMs (see Sec- 
tion |5. 1[ >. Therefore, the best-fit parameters of our model ap- 
pear tobe consistent with other observables. 

The first step to test if this simple model is realistic is to 
check its predictions for RM outside the VLA/CGPS cover- 
age. We do so first by comparing the RM prediction in a strip 
from the south to the north Galact ic pole in the same longi 
tude range. Beyond b= ±30°, the |Taylor et al. 



can be used to test the prediction of our model. In Figure [9] 
we plot the predicted RM as a fun ction of latitude to wards 
Z=108.5° from b=-9Q° to +90°. The jTaylor et "aTJ ( f2009] l RMs 
are binned every 2° and the error bars represent the standard 
deviation of RM within each bin. Our model can reproduce 
the RM-vs-b trend in the observed longitude range from the 
south to the north Galactic pole. We note that the disagree- 
ment between NVSS RMs and the RM prediction at b=+5Q° 
is due to the edge of radio continuum Loop III, w hich could 
contri bute at least +18 rad m" 2 to the RM of EGSs ( |Spoelstra 
|1972| . We did not include a ~ 0.3 /xG vertical magnetic field 
component below the Galactic plane ( |Mao et aL]|2010[ ). We 
note, however, that a constant and small vertical field would 
only shift the entire curve up or down slightly and therefore 
would not affect the overall shape of RM against the b trend. 

Next, we compare the predicted RM from our model with 
the observed RMs in the entire second quadrant from the 
south to the north Galactic pole. This is illustrated in Fig- 
ure|4] our model can reproduce large and negative RMs below 

6 We take notice that since the existing halo field models are global fits to 
the entire high latitude RM sky. rather than a local fit to the RM distribution 
in the region 100°</<117°, 15°<|2>|<30°. Therefore, the quality of fit of 
these global models could be degraded within our specific region of interest 
when attempting to reproduce the RM distribution elsewhere. 



the Galactic plane (Region A). Above the Galactic plane, be- 
sides discrete regions where RMs are affected by foreground 
HII regions and the boundary of the radio continuum Loop III 
(at b ~ +50°), the predicted RM distribution matches the ob- 
served EGS RMs well. The match between the NVSS RMs 
and our model prediction can be further improved by varying 
the Galacto-centric radius at which this non-zero azimuthal 
field exists, specifically by moving the region of non-zero 
magnetic field to larger Galacto-centric radii at larger Galac- 
tic longitudes. This suggests that our halo field model that can 
reproduce the RM-vs-b trend towards 100°</<117° might 
belong to a logarithmic magnetic spiral arm in the Galactic 
halo. In the longitude range 100°<Z<117°, the best fit param- 
eters place a strong coherent magnetic field between Galacto- 
centric radii 8.8 and 10.3 kpc, roughly coinciding with the re- 
gion between the local and Perseus arm (at a Galacto-centric 
radius of ~ 9.95 kpc at Galactic longitude of 134°). This 
could be similar to the magnetic arm phenomenon seen in ex- 
ternal galaxies such as M5 1 and NGC 6946, in which coher- 
ent magnetic fields are obser ved to be s tronger in inter-arm 
regions than within the arms (Beck 2007). 

5.5. Summary and Discussion 



As shown in Sections [5 . 3 1 and pA\ existing exponential and 
double-toroid Galactic halo field models cannot reproduce 
the RMs in our surveyed longitude range 100°</<117° at 
|fe|>15°. A simple double-toroid halo field also does not con- 
form to theory and does not provide a good match to observa- 
tions of the Galactic halo magnetic field. In Section |5~4} we 
show that the presence of a coherent magnetic field within 
radial range 8.8-10.3 kpc at a distance 0.8-2 kpc from the 
Galactic plane can reproduce the RM-vs-b trend within the 
CGPS/VLA survey region. Our model places a strong co- 
herent field within interarm region towards the outer Galaxy. 
Moreover, we find that shifting the radial range of non-zero 
magnetic field strength to larger radii at larger longitude can 
provide a good match to the RM distribution in the entire 
second Galactic quadrant. This is suggestive of a spiral-like 
G alactic halo magnetic fiel d. 

|Beuermann et al.| {1985} used the 408 MHz |Haslam et aL] 
( 1982) map and decomposed the Galactic synchrotron emis- 
sion into a thin disk and a thick disk. The authors estimated 
that the Milky Way synchrotron thick disk has a scale height 
of 1.5 kpc at the location of the Sun and increases to 2.7 kpc 
beyond a Galacto-centric radius of 12 kpc. The synchrotron 
thick disk has comparable height to our halo magnetic field, 
supporting the presence of magnetic fields at distances 1-2 
kpc from the Galactic mid-plane. Moreover, this thick disk 
exhibits spiral structure, and the magnetic field in the thick 
disk is modeled to align moderately with the spiral arms. We 
suggest that RMs in our VLA/CGPS latitude extension region 
could be probing the spiral magnetic field in the thick syn- 
chrotron disk. 

Further insights into this proposed spiral-like thick disk 
field can be obtained by studying magnetic fields in Milky 
Way-ty pe f ace-on galaxies such a s M51. [Berkhuijse n et al.] 
( fl997| ) and |Fletcher et aT] (|20TT) modeled M 51's disk and 
halo field simultaneously, and found spiral magnetic fields (of 
different axial symmetry) in both its disk and its halo. The 
halo field has similar radial extent to the disk field with a to- 
tal Faraday depth of ~ 100 rad m" 2 . This is in agreement 
with our fitted Milky Way halo field: the halo field remains 
substantial beyond the Solar radius, with |RM| up to 100 rad 
m" 2 . In M51, the ratio of the total RM through the disk to 
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that through the halo is ~ 2, which is similar to the observed 
value for the Milky Way. Moreover, similar to what we found 
in the Milky Way towards the outer Galaxy, strong coherent 
halo magnetic field has been found in ga lactocentric radii tha t 
correspond to interarm regions in M5 1 ( [Fletcher et al.|201 l| l. 
Other external galaxies ar e likely to possess spiral-like halo 
fields: |Braun et al.| (2010) modeled the GHz polarized emis- 
sion from a sample or SINGS galaxy with axisymmetric thick 
disk spiral fields with quadrupole topology from the near-side 
out to ~ 30% of the disk radius of these galaxies. This is 
further evidence that the Milky Way could possess magnetic 
fields in the halo at ~ few kpc from the mid-plane that trace 
out spiral arms. 

We propose that modeling of the Milky Way's halo mag- 
netic field might require similar approaches to those used to 
model the magnetic field of the Galactic disk (e.g., Brown 
|et al.||2007]l . Sight lines towards mid- to high- Galactic lati- 
tudes could be divided into spiral arm / inter-arm regions as 
well as plane parallel slabs (in the z direction). Each region 
would have a therm al electron dens i ty de scribed by the WIM 
thick disk model of Gaensler et al. (2008). The strength and 
direction of magnetic fields within these regions can then be 
derived by performing a fit to the observed EGS RMs em- 
pirically. Carrying out such a fitting procedure to the en- 
tire Galactic halo is outside the scope of this paper, but we 
suggest that this will yield a complete picture of the Galactic 
halo magnetic field. Locations of magnetic field reversals as 
well as regions with enhanced magnetic field strength in the 
halo can be identified, as well as the overall axial and vertical 
field symmetry. A better characterization of the Milky Way 
halo magnetic field would be extremely useful to constrain 



its origin. The existing EGS RMs from Taylor et al. (2009) 
combined with those from future surveys, such as the Polar- 
ization Sky Survey of the Universe's Magnetism (POSSUM) 
to be conducted with the Australian Square Kilometre Array 
Pathfinder (ASKAP), will provide an evenly sampled all sky 
RM grid for such studies. 



6. CONCLUSION AND FUTURE WORK 

In this paper, we set out to characterize the symmetry of 
the disk and halo magnetic field of the Milky Way. Build- 
ing on the CGPS latitude extension region, we have con- 
ducted a VLA survey of EGS RMs in the longitude range 



100°</<117 o ,-30 o <Z?<-3 o and+18°<fc< +30°, from which 
we extracted 341 reliable RMs. We have shown that the disk 
field in this region is symmetric about the Galactic mid-plane, 
consistent with the lack of warp of the ionized disk towards 
the outer Galaxy. 

We have shown that towards mid Galactic latitudes in our 
surveyed area, RMs of EGSs are likely produced by a global 
halo magnetic field combined with thermal electrons in the 
thick WIM disk. We demonstrate the unsatisfactory predic- 
tions for the Galactic RM of existing Milky Way halo mag- 
netic field models in the literature. Finally, we have con- 
structed a simple halo field model with reasonable parameters 
that can reproduce the RM distribution in the longitude range 
100°</<117° from b=-90° to +90° reasonably well. This 
model can be extended to explain the EGS RM distribution 
in the entire second Galactic quadrant. We propose future 
modeling of the Galactic halo magnetic field might require 
techniques similar to those used to model the disk magnetic 
field. In future work, we will make predictions of the syn- 
chrotron e miss ion using the simple halo field model proposed 
in Section 5.4 and compare it to existing Galactic synchrotron 
emission survey data. We also plane to extend our simple 
large-scale spiral halo field model to fit to EGS/pulsar RMs at 
mid- to high- Galactic latitudes. 
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FIG. 1.— Comparison between NVSS RMs derived by |Taylor et aT]j2009} and our RMs in 100°</<117°, |£>|<30°. The solid line of slope 1 indicates where 
NVSS RMs and our VLA RMs agree with each other. 
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FIG. 2. — RM distribution in the region 100°</<117°, |fc|<30° overlaid on an Ha intensity map I Finkbeiner 2003). The color scale is in units ofRayleighs. 
Positive (negative) RMs are denoted by filled (open) circles with diameters proportional to |RM|. Sources with KMs consistent with zero at Icr are denoted by 
asterisks. 
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FIG. 3 . — Top: The variation of RM as a function of Galactic latitude for EGSs in the Galactic longitude range 100° </< 1 17° . The RM data are binned such 
that there are 40 RMs within each bin. We have plotted the median RM and the standard error of the mean within each bin. Bottom: The variation of RM as a 
function of the absolute value of the Galactic latitude for EGSs in the Galactic longitude range 100°</<117°. This plot is the same as in the top panel except 
that the data have been folded across b=0°. Black (red) symbols represent the RM-vs-b behavior above (below) the Galactic plane. 




Galactic Longitude 

FIG. 4. — The RM predicted from our simple halo field model towards the entire second Galactic quadrant viewed as a 2D distribution: the col or scheme is 
chosen such that red (blue) represents positive (negative) RMs, while white represents RMs close to rad m~ 2 . RMs of extragalactic sources from |Taylor et al.| 
(2009) are over-plotted as filled circles on the same color scheme. Three large regions where there is discrepancy between the model fit and the KM data are 
labeled: Loop III and HII regions Sivan 3 and Sh 2-220. The Milky Way disk magnetic field dominates within 15° of the Galactic plane and therefore the region 
is blanked. 
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FIG. 5, — Toroidal magneti c field strength in a vertical cross-section of the Milky Way using the Jansson et al. 1 2009) best fit double torus model (top panel) 
and the Ruiz-Granados et al. (20101 model (bottom panel). 
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FIG. 6. — The trend of RM as a function of Galactic latitude towards /=108.5°. RM predictions from 6 different halo field models in the literature have been 
plotted. We assume that the halo magnetic field dominates at a height of 540 pc, which provides an upper limit for |RM| . (Models with larger transition heights 
have fewer thermal electrons and hence smaller |RM|). The VLA and CGPS RMs are binned every 3°: diamonds represent the median RM within each bin 
and the error bars denote the weighted standard deviation within each bin. For |fo|<15° (shaded region), the disk field dominates and therefore is not shown for 
clarity. 
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FIG. 7. — The trend of RM as a function of Galactic latitude towards Z=108.5°. The RM prediction from our halo field model in Section |5!4] (solid line) is 
over-plotted on the binned VLA and CGPS RMs denoted by diamonds. The error bars are the weighted standard deviation within each bin. For |to|<15° (shaded 
region), the disk field dominates. 
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GLON (degrees) 

Galactic Longitude 

FIG. 8. — Left Panel: The RM predicted from our simple halo field model shown as a 2D distribution: the color scheme is chosen such that red (blue) represents 
positive (negative) RMs, while white represents RMs close to rad irT 2 . Measured RMs of EGSs are overplotted as filled circles on the same color scale. 
Color differences between filled circles and the background shows the discrepancy between the model prediction and the observed RMs. We note that the 
uncertainties associated with individual EGS RMs are not represented in this image. The color difference between the observed and predicted RM in latitude 
range +15°<fc<+20° is likely due to the larger measurement errors of the CGPS RMs compared to the VLA RMs. Right Panel: The di stribution of the re sidual 
RM (the difference between the predicted RMs from our halo field model and the measured RMs) overlaid on an Ho intensity map iFinkbeiner 2003). The 
color scale is in units of Rayleighs. Positive (negative) residual RMs are denoted by filled (open) circles with diameters proportional to the amount of residual. 
Asterisks represent sight lines with residual RMs consistent with zero at lcr. 
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FIG. 9. — The trend of RM as a function of Galactic latitud e towards /= 108.5°. The predicted RM against b trend from our halo field model in Section 
(solid line) is over-plotted on the binned NVSS RMs I Taylor et al. 2009) denoted by diamonds. Error bars represent the weighted standard deviation of R 
within each bin. Since the Milky Way disk field dominates within 15" ot the Galactic plane, RMs in this region are not plotted for clarity. 
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Rotation Measures of Extragalactic Sources from VLA observations in 
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+ 38:13:48.00 


106.89 


-22.27 


-55±6 


18:08:14.40 


+ 75:49:12.00 


106.99 


+28.88 


-66± 3 


23:24:02.40 


+ 44:55:12.00 


107.08 


-15.24 


-23± 9 


23:31:07.20 


+ 41:16:48.00 


107.08 


-19.12 


-36±4 


23:01:55.20 


+ 53:47:00.00 


107.09 


-5.67 


-49±12 


23:39:28.80 


+ 36:20:12.00 


107.09 


-24.31 


-86±10 


23:14:31.20 


+ 49:17:12.00 


107.13 


-10.59 


-27± 7 


23:40:45.60 


+ 35:53:60.00 


107.21 


-24.82 


-128± 9 


23:34:28.80 


+ 39:46:36.00 


107.22 


-20.73 


-51±10 


23:21:24.00 


+ 46:34:12.00 


107.23 


-13.53 


- Bill 


18:30:09.60 


+ 76:06:36.00 


107.30 


+27.55 


-66± 8 


23:20:57.60 


+ 47:05:48.00 


107.34 


-13.02 


- 11± 6 


23:06:57.60 


+ 52:45:60.00 


107.36 


-6.93 


- 37±20 


23:36:52.80 


+ 38:50:00.00 


107.38 


-21.77 


-60±17 


23:29:43.20 


+ 43:04:36.00 


107.44 


-17.33 


-34±9 


23:05:28.80 


+ 53:30:24.00 


107.46 


-6.15 


-113±13 


23:18:26.40 


+ 48:37:48.00 


107.49 


-11.44 


-7± 6 


23:43:50.40 


+ 34:49:24.00 


107.54 


-26.01 


-63± 8 


23:06:16.80 


+ 53:31:24.00 


107.57 


-6.19 


-75±11 


23:41:12.00 


+ 36:43:36.00 


107.58 


-24.04 


-43±9 


23:38:52.80 


+ 38:23:60.00 


107.64 


-22.32 


-58±4 



TABLE 1 — Continued 



t? A/T9nnnvvimsA 
Krt^j zuuu ) ^nms ) 




1 K ) 


\ ) 


Kivi ^rdti rn ) 


20:17:14.40 


+ 74:40:00.00 


107.69 


+20.72 


-27± 5 


17:55:55.20 


+ 76:27:00.00 


107.79 


+29.56 


-25±4 


23:14:19.20 


+ 51:10:12.00 


107.81 


-8.83 


-101±13 


23:48:50.40 


+ 32:08:12.00 


107.81 


-28.88 


-70±3 


23:08:45.60 


+ 53:21:00.00 


107.85 


-6.49 


-80±14 


23:38:09.60 


+ 39:33:48.00 


107.87 


-21.16 


-54±7 


23:17:12.00 


+ 50:23:12.00 


107.95 


-9.72 


- 10±11 


23:26:57.60 


+ 46:19:24.00 


108.07 


-14.09 


+ 12±10 


23:11:14.40 


+ 53:10:12.00 


108.12 


-6.79 


-78±10 


19:34:36.00 


+ 76:17:12.00 


108.20 


+23.80 


-32±14 


18:45:02.40 


+ 76:52:36.00 


108.21 


+26.73 


-84±2 


23:15:07.20 


+ 52:06:24.00 


108.27 


-8.01 


-65± 7 


19:04:26.40 


+ 76:47:12.00 


108.27 


+25.62 


-58±11 


23:33:21.60 


+ 43:46:00.00 


108.33 


-16.87 


-35± 6 


23:46:33.60 


+ 35:28:00.00 


108.33 


-25.56 


-65±4 


23:32:04.80 


+ 44:30:12.00 


108.34 


-16.11 


- 14± 4 


23:23:07.20 


+ 49:05:48.00 


108.40 


-11.27 


-61±30 


23:41:14.40 


+ 39:18:24.00 


108.40 


-21.57 


-62±12 


23:25:45.60 


+ 48:06:12.00 


108.47 


-12.34 


- 44± 1 1 


23:16:21.60 


+ 52:19:24.00 


108.53 


-7.88 


-81±9 


23:35:48.00 


+ 43:06:48.00 


108.55 


-17.65 


-33±11 


23:51:43.20 


+ 32:21:48.00 


108.55 


-28.83 


-68±10 


23:31:48.00 


+ 45:22:36.00 


108.58 


-15.26 


- 12± 9 


18:43:14.40 


+ 77:15:24.00 


108.62 


+26.86 


-75±9 


23:37:24.00 


+ 42:28:12.00 


108.65 


-18.34 


-43± 8 


19:32:04.80 


+ 76:46:36.00 


108.66 


+24.08 


-41±10 


23:27:26.40 


+ 48:02:12.00 


108.73 


-12.49 


-22± 4 


23:25:31.20 


+ 49:01:48.00 


108.74 


-11.47 


-47±19 


20:38:07.20 


+ 75:05:48.00 


108.80 


+19.75 


-18±7 


23:48:21.60 


+ 35:54:12.00 


108.85 


-25.23 


-82±6 


23:53:52.80 


+ 32:09:00.00 


108.99 


-29.15 


-54±9 


23:25:38.40 


+ 50:03:12.00 


109.12 


-10.50 


-73±11 


23:39:19.20 


+ 42:57:00.00 


109.15 


-17.98 


- 65±20 


23:55:43.20 


+ 31:25:48.00 


109.22 


-29.95 


-98±19 


23:52:38.40 


+ 34:10:12.00 


109.28 


-27.14 


-48±23 


20:22:36.00 


+ 76:11:24.00 


109.29 


+21.13 


-34± 3 


23:42:48.00 


+ 41:16:12.00 


109.31 


-19.77 


-512±32 


23:28:50.40 


+ 49:19:36.00 


109.37 


-11.36 


-76± 8 


23:49:43.20 


+ 36:59:12.00 


109.45 


-24.26 


-112±14 


23:53:19.20 


+ 34:17:36.00 


109.48 


-27.06 


-90±15 


18:07:40.80 


+ 78:04:48.00 


109.56 


+28.75 


-64±12 


20:13:57.60 


+ 76:46:36.00 


109.60 


+21.85 


-35±10 


20:56:43.20 


+ 75:14:12.00 


109.64 


+18.85 


-35± 2 


23:26:04.80 


+ 51:40:48.00 


109.73 


-8.99 


-155±12 


23:49:21.60 


+ 38:49:36.00 


109.89 


-22.47 


-103±19 


23:56:48.00 


+ 33:04:36.00 


109.93 


-28.42 


-84± 8 


23:40:07.20 


+ 45:34:36.00 


110.08 


-15.52 


-62±10 


20:12:48.00 


+ 77:19:60.00 


110.09 


+22.15 


-44±4 


23:42:12.00 


+ 44:41:00.00 


110.19 


-16.47 


-66±17 


23:26:31.20 


+ 52:57:24.00 


110.22 


-7.80 


-123±22 


23:39:09.60 


+ 46:41:12.00 


110.24 


-14.40 


-45±16 


20:45:40.80 


+ 76:25:00.00 


110.25 


+20.08 


- 12± 2 


23:59:50.40 


+ 31:40:48.00 


110.27 


-29.93 


-96±15 


23:34:43.20 


+ 49:15:36.00 


110.28 


-11.72 


-35± 7 


23:30:07.20 


+ 51:54:48.00 


110.40 


-8.97 


-514±25 


23:41:36.00 


+ 45:49:24.00 


110.41 


-15.34 


-43±12 


23:41:36.00 


+ 46:03:12.00 


110.48 


-15.11 


-50± 4 


23:48:45.60 


+ 41:27:48.00 


110.50 


-19.90 


-61±16 


21:18:14.40 


+ 75:12:12.00 


110.51 


+17.76 


+ 7± 1 


21:27:48.00 


+ 74:50:48.00 


110.67 


+17.06 


- 12±12 


23:26:24.00 


+ 54:49:36.00 


110.80 


-6.03 


-98±12 


23:36:26.40 


+ 50:08:24.00 


110.82 


-10.95 


- 60±20 


23:44:36.00 


+ 46:04:12.00 


111.00 


-15.25 


- 46±30 


21:22:57.60 


+ 75:37:24.00 


111.05 


+17.82 


+ 25± 3 


23:46:24.00 


+ 45:19:12.00 


111.11 


-16.05 


-438±16 


19:43:38.40 


+ 78:58:36.00 


111.15 


+24.15 


- 16± 1 


20:12:40.80 


+ 78:28:12.00 


111.22 


+22.65 


+ 8±7 


23:40:52.80 


+ 49:36:24.00 


111.36 


-11.68 


-482±27 


23:35:50.40 


+ 52:15:12.00 


111.36 


-8.91 


+336±11 


23:50:38.40 


+ 43:40:12.00 


111.44 


-17.84 


-60±10 


20:08:52.80 


+ 78:55:00.00 


111.57 


+23.00 


+ 10±12 


21:02:36.00 


+ 77:15:12.00 


111.57 


+19.78 


-14±13 


23:49:21.60 


+ 45:35:36.00 


111.70 


-15.92 


- 50±24 


18:53:09.60 


+ 80:01:00.00 


111.77 


+26.61 


- 15±29 


23:44:12.00 


+ 49:07:48.00 


111.77 


-12.28 


-41±21 


00:01:55.20 


+ 36:22:48.00 


111.92 


-25.44 


-104±13 


23:47:07.20 


+ 48:01:36.00 


111.95 


-13.47 


-33±17 
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TABLE 1 — Continued 



Krt^j zuuu ) (nms ) 




1 K ) 


\ ) 


R1VT CrctA m _ 2\ 
K1V1 ^rdti m ) 


21:19:57.60 


+ 76:57:12.00 


111.98 


+18.85 


+ 2±9 


23:50:16.80 


+ 46:06:36.00 


112.00 


-15.46 


-46±18 


17:43:40.80 


+ 80:04:00.00 


112.00 


+29.60 


-42± 4 


19:34:19.20 


+ 79:56:00.00 


112.03 


+24.84 


+ 14± 3 


23:49:04.80 


+ 47:21:12.00 


112.11 


-14.20 


- 16±15 


23:38:52.80 


+ 53:37:60.00 


112.19 


-7.73 


-108±20 


19:12:24.00 


+ 80:26:36.00 


112.34 


+25.86 


-20±5 


23:37:45.60 


+ 55:51:12.00 


112.67 


-5.55 


-109±32 


23:53:28.80 


+ 46:51:24.00 


112.73 


-14.87 


-52±22 


00:01:28.80 


+ 41:04:24.00 


112.89 


-20.82 


-74±14 


20:58:43.20 


+ 79:05:00.00 


113.01 


+21.03 


+ Hill 


00:00:38.40 


+ 43:57:48.00 


113.36 


-17.97 


-70± 8 


23:52:12.00 


+ 50:26:24.00 


113.37 


-11.33 


- 1±22 


21:57:12.00 


+ 76:46:00.00 


113.43 


+17.24 


- 19± 5 


18:11:48.00 


+ 81:30:24.00 


113.43 


+28.27 


+ 26±4 


21:54:07.20 


+ 77:04:48.00 


113.49 


+17.57 


+ 10± 4 


23:51:00.00 


+ 51:49:24.00 


113.52 


-9.94 


-55±21 


00:05:07.20 


+ 40:57:60.00 


113.58 


-21.08 


-72±23 


23:47:28.80 


+ 54:32:12.00 


113.67 


-7.17 


-38±21 


00:11:52.80 


+ 34:16:36.00 


113.69 


-27.88 


-73±13 


18:45:33.60 


+ 81:50:24.00 


113.76 


+27.04 


+ 15± 1 


00:03:19.20 


+ 43:34:36.00 


113.77 


-18.44 


-51± 8 


17:51:00.00 


+ 81:47:60.00 


113.87 


+28.98 


- 17± 4 


00:12:40.80 


+ 35:39:12.00 


114.14 


-26.55 


-76±33 


00:14:45.60 


+ 33:40:36.00 


114.24 


-28.59 


-65±18 


00:04:55.20 


+ 44:27:48.00 


114.25 


-17.63 


-59±12 


21:09:16.80 


+ 80:20:24.00 


114.40 


+21.42 


+ 21±14 


00:06:14.40 


+ 44:04:24.00 


114.42 


-18.04 


-41± 8 


23:51:40.80 


+ 55:32:60.00 


114.49 


-6.35 


-85±17 


00:15:14.40 


+ 35:03:48.00 


114.60 


-27.23 


-75± 6 


21:15:50.40 


+ 80:29:00.00 


114.69 


+21.29 


+ 31± 6 


19:38:36.00 


+ 82:33:36.00 


114.87 


+25.38 


+ 14±28 


00:15:31.20 


+ 36:12:36.00 


114.87 


-26.11 


-82±12 


23:57:28.80 


+ 53:31:24.00 


114.87 


-8.50 


-53±25 


23:59:43.20 


+ 52:37:36.00 


115.02 


-9.45 


+ 10±10 


00:12:55.20 


+ 40:32:24.00 


115.07 


-21.73 


- 80±12 


23:56:38.40 


+ 55:18:24.00 


115.13 


-6.73 


-55±15 


00:11:07.20 


+ 43:15:48.00 


115.19 


-19.00 


- 15±10 


23:57:00.00 


+ 55:41:48.00 


115.26 


-6.37 


-69±5 


19:45:28.80 


+ 82:57:00.00 


115.35 


+25.28 


+ 8±6 


00:18:52.80 


+ 34:38:12.00 


115.37 


-27.75 


-51±11 


23:55:33.60 


+ 57:28:60.00 


115.45 


-4.59 


-66±11 


21:10:07.20 


+ 81:36:36.00 


115.49 


+22.16 


-5±10 


00:12:14.40 


+ 43:58:36.00 


115.52 


-18.32 


-49±10 


21:25:04.80 


+ 81:15:36.00 


115.57 


+21.51 


+ 22±36 


17:09:16.80 


+ 83:00:12.00 


115.65 


+29.95 


-48±5 


00:06:07.20 


+ 50:49:36.00 


115.66 


-11.40 


+ 37± 4 


22:17:19.20 


+ 78:58:24.00 


115.72 


+18.32 


-21±8 


17:05:07.20 


+ 83:16:12.00 


115.98 


+29.97 


-40± 3 


22:15:57.60 


+ 79:29:12.00 


115.99 


+18.78 


-26± 6 


00:17:33.60 


+ 40:47:12.00 


116.04 


-21.63 


-74±5 


00:10:40.80 


+ 48:47:48.00 


116.05 


-13.52 


-9±9 


00:12:40.80 


+ 47:04:48.00 


116.11 


-15.29 


- 18± 8 


00:11:00.00 


+ 49:04:24.00 


116.15 


-13.26 


-8±4 


00:05:50.40 


+ 54:00:00.00 


116.19 


-8.28 


-49±17 


22:18:24.00 


+ 79:41:24.00 


116.21 


+18.88 


-8±3 


00:05:16.80 


+ 55:17:12.00 


116.34 


-6.99 


-53±14 


00:07:24.00 


+ 53:32:60.00 


116.34 


-8.77 


-2±5 


00:11:52.80 


+ 49:31:00.00 


116.37 


-12.84 


+ 18±12 


00:04:28.80 


+ 56:14:12.00 


116.40 


-6.04 


-56±7 


21:27:50.40 


+ 82:13:24.00 


116.42 


+22.07 


+ 16±5 


00:22:12.00 


+ 36:57:60.00 


116.46 


-25.55 


-104± 8 


00:24:43.20 


+ 32:51:12.00 


116.48 


-29.68 


-46±19 


18:31:43.20 


+ 84:16:60.00 


116.49 


+27.47 


+ 4± 7 


00:11:36.00 


+ 50:41:60.00 


116.51 


-11.68 


+ 18± 6 


21:09:38.40 


+ 82:55:48.00 


116.60 


+22.97 


-8±10 


00:15:16.80 


+ 47:36:48.00 


116.64 


-14.82 


-37± 5 


00:25:19.20 


+ 33:00:00.00 


116.64 


-29.54 


-50± 8 


00:16:43.20 


+ 46:08:48.00 


116.68 


-16.31 


-20± 8 


20:23:33.60 


+ 83:54:12.00 


116.77 


+24.61 


-15±18 


00:11:14.40 


+ 53:39:12.00 


116.92 


-8.75 


-21±6 


00:15:40.80 


+ 49:07:48.00 


116.93 


-13.33 


-32±3 


00:19:14.40 


+ 44:49:36.00 


116.95 


-17.67 


-33± 6 



TABLE 2 

Double-Torus Halo Magnetic Field Parameters 



Reference 


Bq (m G) 


( k P c ) 


Zo (kpc) 


zf (kpc)for| z |<z« 


zf (kpc) for |z|>*» 


Model type 


Sunetal.||2008) 


10 


4 


1.5 


0.2 


0.4 


Qualitative 


Jansson et al. IZ009I 


4.9 


18 


1.4 


0.12 


8.5 


Fitted" 


Sun & KeicHf(2UlUf 


2 

4 for z>0 
2 for z<0 


4 


1.5 


0.2 


4 


Qualitative 


Pshirkovetal.||2011f 


6 


1.3 


0.25 


0.4 


Benchmark b 



These parameters were found by fitting to EGS RMs and diffuse polarized synchrotron emission. 
These parameters were chosen based on model fits to EGS RMs. 



